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Starting from Stratton-Panofsky-Phillips-Jefimenko equations for the electric and magnetic fields generated 
by completely arbitrary charge and current density distributions at rest, we derive far-zone approximations for 
the fields, containing all components, dominant as well as sub-dominant. Using these approximate formulas, we 
derive general formulas for the total electromagnetic linear momentum and angular momentum, valid at large 
distances from arbitrary, non-moving charge and current sources. 



I. INTRODUCTION 

Standard classical electrodynamics offers different meth- 
ods for calculating electromagnetic fields and physical observ- 
ables derived from them. Usually the methods amount to solv- 
ing the Maxwell equations in one form or another, or to find 
the potentials, choose a gauge and then perform the calcula- 
tions. 

In recent years, the use of alternative methods, based on in- 
tegral equations that relate the fields directly to their charge 
and current sources, have gained attraction; see Ref. 1 and ref- 
erences cited therein. Here we show how these methods can 
be used to directly derive general expressions for the linear 
and angular momentum of the electromagnetic field that ap- 
proximate these quantities at large distances from the sources. 
Among other things, these expressions show that both linear 
momentum and angular momentum can be used to transfer 
information wirelessly in free space. 

In Chapter VIII of his 1941 textbook, Stratton^ calcu- 
lates, without the intervention of potentials, temporal Fourier 
transform expressions for the retarded electric and magnetic 
fields, E(f,x) and B(r,x), respectively, generated by arbi- 
trary distributions of charge and current densities at rest rel- 
ative to the observer. In Chapter 14 of the second edi- 
tion of their textbook on electrodynamics, published in 1962, 
Panofsky and Phillips^ present a variant form of these ex- 
pressions, and give them also in ordinary space-time coor- 
dinates. Four years later, Jefimenko published his electro- 
dynamics textbook'* where the Panofsky and Phillips expres- 
sions for the retarded E and B fields were given in Chap- 
ter 15. These expressions are sometimes referred to as the 
Jefimenko equations, ^'^ but should, perhaps, rather be called 
the Stratton-Panofsky-Phillips-Jefimenko (SPPJ) equations, a 
name we will use throughout. 

Introducing the notation x for the observer's coordinate, x' 



for the source coordinate, and 



x-x' 



(1) 



for the retarded time' relative to the source point where t is the 
observer's time and c is the speed of light, the SPPJ equations 
for the retarded electric and magnetic fields can, in obvious 
notation, be written 
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FIG. 1. The sources are located at points x' inside a volume V' lo- 
cated around xq such that sup |x' — xq | <C inf |x — x' | , where x is the 
observation point. The unit vector k is directed along x — x'. n is 
along X — X() and orthogonal to the surface S{x()) with its origin at xq. 



the SPPJ equations (2) can be cast into the more compact and 
symmetric form 



(4a) 



so that Eq. (2a) can be written in the alternative form 
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which more clearly exhibits the relation between the various 
components of the retarded E field. 

In Section II, far-zone formulas for the E and B fields, based 
on the SPPJ equations, are introduced. These equations are 
used in Section III and IV to derive far-zone expressions for 
the electromagnetic linear momentum and angular momen- 
tum, respectively. In Section V, a summary and conclusions 
are made. 



II. THE FIELDS AT LARGE DISTANCES FROM THE 
SOURCE 

Following Panofsky and Phillips ^ we introduce the tempo- 
ral Fourier component representations of Eq. (5) for E(f,x) 



and ofEq. (2b) for B(f ,x). 
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respectively. 

As illustrated in Fig. 1, the observation point x is as- 
sumed to be located far away from the sources, which in 
turn are assumed to be localized near a point xo inside 
a volume V' that has such a limited spatial extent that 
sup |x' — xol ^ inf |x — x'|, and the integration surface S, cen- 
tered on Xq, has a large enough radius |x — xo| ^ sup |x' — xo|. 
Then one can make the usual far-zone approximation 
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so that Eq. (5) for the electric field and Eq. (2b) for the mag- 
netic field can, in complex notation, be approximated by 
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as obtained from inverse Fourier transforming of Eqs. (6) with 
Eq. (7) inserted. 

In many cases, when one wants to calculate electromagnetic 
observables in the far zone, the approximate Eqs. (10) are ac- 
curate enough. At the same time they are easier to evaluate, 
and give results which are physically more lucid than the SPPJ 
equations (2). Below, we follow this scheme to derive general, 
untruncated expressions for the linear and angular momentum, 
valid far away from the source volume Y' . 



III. LINEAR MOMENTUM 



The linear momentum conservation law^ 
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describes the balance between the time rate of change of the 
mechanical linear momentum p™'^'^^^ ^ mechanical force, 
the time rate of change of the electromagnetic field linear mo- 
mentum p'^'^''*, and the flow of linear momentum across the sur- 
face S described by the linear momentum flux tensor density 
T (the negative of Maxwell's stress tensor) of the electromag- 
netic field. For example, this law describes how the transla- 
tional motion of charges (e.g., a transmitting antenna current) 
at one point in space induces a translational motion, via force 
action, of remote charges {e.g., a receiving antenna current). 
This is one of several physical mechanisms by which informa- 
tion can be transferred electromagnetically through free space 
and the primary physical basis for current radio astronomy and 
other radio-based science, as well as today's wireless commu- 
nications technology. 

The linear momentum of the electromagnetic field is de- 
fined as^'^ 



, field 



where 
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is the electromagnetic field linear momentum density and S 
the Poynting vector. 

Applying the paraxial approximation^ 
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of the wave vector direction, and moving the constant unit 
vector n, from xq to the observation point x, outside each of 



the six constituent integrals in in the following way 

^dVp(f',x')k« (^dVp(f',x'))n = ^(f')n (15a) 



dYj(f',x') xkw / dVj(f',x') xn = I(f')xn 



(15b) 



dVj(f',x') xk: 



dV[j(f',x')-k]k: 



dVj(f',x') xn = t(f') xn 



dVj(f',x'))-fl 



= (I-n)n=/„(Ofi 



dV[j(f',x') xk] xk! 



dVj(f',xO) xn 



(15c) 



(15d) 



X n 



(Ixn) xn = /„(f')n-I(r') 



(15e) 



dV[j(f',x') xk] xk; 



dVj(f',x')) xn 



X n 



(ixn) xn = /„(f')n-i(r') 

(15f) 



one obtains, after some vector algebraic simplifications, the 
following general, untruncated expression for the cycle aver- 
aged linear momentum density 
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This approximate expression is valid far away from the source 
volume y'. 

At very large distances r = |x — xo| from the source vol- 
ume y', we see that the linear momentum density, and conse- 
quently also the Poynting vector, is accurately represented by 
the first term, which is radial (along n) and falls off as 1 / r^, 
while the successive non-radial terms fall off as Xjr' or as 
1 / r^, respectively. However, at finite distances the Poynting 
vector is not radial but has transverse components. Of course, 
these transverse components, which make the Poynting vector 
spiral around n, are difficult to observe because of their small- 
ness. Furthermore, since the non-radial terms fall off faster 
than 1 / r-, the spiraling of the Poynting vector will diminish 
with distance from the source. 
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IV. ANGULAR MOMENTUM 



The angular momentum conservation law 
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describes the balance between the time rate of change of the 
mechanical angular momentum J™'^'^'*, i.e., mechanical torque, 
the time rate of change of the electromagnetic field angular 
momentum J*''^''^, and the flow of angular momentum across 
the surface S described by the angular momentum flux ten- 
sor K = (x — xo) X T, all about the point xq. This law de- 
scribes how the rotational (spin, orbital) motion of charges at 
one point in space induces a rotational motion, via torque ac- 
tion, of charges at other points in space. This is an additional 
physical mechanism by which information can be transferred 
electromagnetically through free space but one that is used 
only partially and sparingly in today's radio-based research 
and wireless communication technology. 

The field angular momentum about a point xo is defined 
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is the electromagnetic angular momentum density. For a sin- 
gle temporal Fourier component in complex notation and a 
beam geometry, Eq. (18) can be written'*^-^' 
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The first term is the electromagnetic orbital angular momen- 
tum (OAM), which describes the vorticity of the EM field, 
and the second term is the electromagnetic spin angular mo- 
mentum (SAM), which describes wave polarization (left or 
right circular). 

Since x — xo = |x — xo|n, the vector product in Eq. (19) 
gives vanishing contributions when operating on the radial 
terms (terms parallel to n) in Eq. (16). As a result, the com- 
plete cycle averaged far-zone expression for a frequency com- 
ponent (0 of the electromagnetic angular momentum density 
generated by arbitrary charge and current sources is simply 
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We see that at very large distances r, the angular momentum 
density falls off as 1/r^, i.e., it has precisely the same behav- 
ior in the far zone as the linear momentum density and can 
therefore also transfer information wirelessly over large dis- 
tances. The only difference is that while the direction of the 
linear momentum (Poynting vector) becomes purely radial at 
infinity, the angular momentum becomes perpendicular to the 
linear momentum, i.e. purely transverse, there. 



V. SUMMARY AND CONCLUSIONS 

We have shown how the SPPJ equations for the retarded 
electric and magnetic fields, generated by arbitrary charge and 
current distributions, can be used to derive general far-zone 
formulas for the concomitant electromagnetic linear and an- 
gular momentum. 

From the far-zone approximations obtained for the linear 
momentum and angular momentum densities, Eq. (16) and 
Eq. (22), respectively, it is easy to see that both these physical 
quantities, to leading order, fall off as 1 /r^ with distance r 
from the source. Consequently, when integrated over a surface 
element sin(0)d0d(p of a large spherical shell, centered on 
the source, they behave as constants. Physically, this means 
that both linear and angular momentum can be carried all the 
way to infinity. While the time rate of change of the linear 
momentum provides the force that causes the charges in the 
EM sensor {e.g., a receiving antenna) perform translational 
(oscillating) motions, the angular momentum gives rise to a 
rotational (spinning and/or orbiting) motion. 

However, while the linear momentum is in the far zone de- 
termined essentially by the dominant E and B far-zone fields 
that fall off as the far-zone angular momentum is not de- 
termined by far-zone fields, but rather by near-zone compo- 
nents of the E field, that fall off as r^^, and the far-zone B 
field. In fact, the far-zone E field does not contribute to the 
far-zone angular momentum at all ! This somewhat surprising 
result is a generalization to arbitrary EM fields of a result, orig- 
inally derived by Abraham already in 1914, for the special 
case of pure dipole fields. To quote from page 916 of Abra- 
ham (using modern notation):''^ 

"Der elektrische Vektor E hat (s. Gl. 6c un- 
ten) eine radiale Komponente, die freilich mit 
wachsendem r wie r^^ abnimmt, wahrend die 
zum Fahrstrahl senkrechten Komponenten von 



E 



abnehmen. Mit wachsender Ent- 



fernung von der Lichtquelle werden somit die 
Lichtwellen transversal, der Impuls wird paral- 
lel dem Fahrstrahl; es konnte darum nach (4) 
scheinen, als ob der Drehimpuls der Wellenzone 
gleich null sei. Doch sieht man, daB dem nicht so 
ist, wenn man die Ordnung der fraglichen GroBen 
bestimmt; wenn auch die longitudinale Kompo- 
nente von E von der Ordnung r^^ ist, so ist doch 
das skalare Produkt von r und E von der Ordnung 
r^', ebenso wie H. Aus (4c) folgt demnach, daB 
die Dichte des Impulsmomentes, ebenso wie die 
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Dichten des Impulses und der Energie, von der 
Ordnung r^^ ist." 

The fact that electromagnetic linear and angular momen- 
tum can be transferred wirelessly over long distances (they 
are both irreversibly lost at infinity), means, e.g., that both of 
them can be used in wireless communications. Today, only 
the linear momentum and, in some cases the two orthogonal 
states of the spin part of the angular momentum (SAM, wave 
polarization), are used in wireless communication protocols, 
including such modern concepts as MIMO.''* However, wave 
polarization techniques do not always work reliably in real- 
world communication settings. A common problem is the de- 
polarization of the radio beam, due to reflections and other 
interactions with the surroundings that can cause a conversion 
of polarization (SAM) into orbital angular momentum (OAM). 
The SAM is thereby diminished or lost, while the total an- 
gular momentum (SAM plus OAM) is still conserved. The 
use of OAM, in addition to SAM (polarization), could there- 
fore be a remedy for this deficiency. Besides, and perhaps 
more importantly, OAM spans a denumerably infinite Hilbert 
space, providing, in principle, an infinite number of orthog- 
onal basis states that can be used as the "letters" of a higher 



dimensional information alphabet. ' By adding OAM en- 
coding to an EM beam, this beam can therefore transfer more 
wireless information per unit time and unit frequency than a 
radio beam carrying only linear momentum (power), or linear 
momentum plus SAM (polarization) as is the case for current 
communication techniques. '^''^ 

The possibility to use OAM encoding of EM beams for ef- 
ficient free-space information transfer has been successfully 
demonstrated in several experiments, ^'^"^^ even at the single- 
photon level. A scheme for the use of OAM at radio fre- 
quencies in the gigahertz range and below, where digital radio 
methods are readily applicable to make the method feasible 
for practical technology applications, has been developed. 
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